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摘要 共 轨 运动 天 体 与 摄 动 天 体 的 半 长 径 相 同 , 处 于 1:1 平 运动 共振 中 . 太阳 系 内 多 个 行星 的 特洛伊 天 体 即 为 处 
于 出 昱 形 轨 道 的 共 轨 运动 天 体 , 其 中 一 些 高 轨道 倾角 特洛伊 天 体 的 轨道 运动 与 来 源 仍 未 被 完全 理解 . 利用 一 个 
新 发 展 的 适用 于 处理 1:1 平 运动 共振 的 摄 动 函 数 展开 方式 , 对 三 维 空间 中 的 共 轨 运动 进行 考察 , 计算 不 同 初 始 轨 
道 根 数 情况 下 共 轨 轨道 的 共振 中 心 、 共 振 宽 度 , 分 析 轨 道 类 型 与 初始 轨道 根 数 的 关系 . 并 将 分 析 方 法 所 得 结果 
与 数值 方法 的 结果 相互 比较 验证 , 得 到 了 广阔 初始 轨道 根 数 空间 内 共 轨 运动 的 全 局 图 景 . 
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1 引言 经 有 10000 颗 左右 的 木星 特洛伊 、32 颗 海王 星 特 洛 
圆 型 限制 性 三 体 问题 存在 五 个 平 动 点 ， 均 位 伊 、9 颗 火星 特洛伊 、1 果 天王星 特 洛 伊 以 及 2 颗 地 
于 两 个 主 天 体 m、m 的 轨道 面 zy 平面 内 . 其 中 两 个 球 特洛伊 被 观测 到 |. 

平 动 点 [4、Ls 分 别 与 m、m/ 构 成 等 边 三 角形 的 构 FEF A AE SIA, AAEE K 
型 , 称 为 三 角 平 动 点 ， 当 系统 的 质量 参数 /, = m// ”与 第 二 主 天 体 相同 ， 因 而 其 运动 被 称 为 “ 共 轨 运 
(m 十 m/ 在 Routh 临 界 质量 (~ 0.03852) 以 下 时 , Las ay”. 常见 的 共 轨 运动 轨道 包括 围绕 某 个 三 角 平 动 
Ls 点 附近 天 体 的 运动 是 线性 稳定 的 趾 , 而 其 非 线性 。 ”在 (4 或 Ls) 附近 的 暗 旦 形 轨 道 、 围 绕 两 个 三 角 平 
稳定 性 ( 仅 除去 = 0.02429、0.01351 两 个 不 稳定 的 动 点 和 一 个 共 线 平 动 点 (Zs) 的 马蹄 形 轨道 以 及 在 
质量 参数 值 ) 则 在 KAM 定 理 (Kolmogorov-Arnold- 第 二 主 天 体 附近 的 准 卫 星 轨道 等 , Hp Aee E 
Moser theorem) 出 现 之 后 得 到 了 证 明 已 色 . 1906 年 ， 道 最 为 常见 . 我 们 在 图 1 中 显示 了 这 几 类 轨道 类 型 . 
天 文学 家 在 木星 轨道 处 观测 到 第 一 颗 在 Ly 点 附近 ”图 中 选用 了 旋转 坐标 系 , 其 旋转 速度 与 主 天 体 相 互 
运动 的 小 天 体 ， 上述 理论 因此 获得 天 文 观测 的 证 。“” 绕 转 的 圆 轨道 角速度 相同 . 图 1 中 的 几 个 例子 , 小 天 
Sk. 后 来 ， 此 类 在 太阳 -行星 的 三 角 平 动 点 附近 运 体 与 两 个 主 天 体 轨 道 共 面 , 而 且 轨 道 的 初始 偏心 率 
动 的 天 体 被 称 为 特洛伊 小 天 体 ( 本 文 以 下 简称 为 都 较 小 . 在 小 天 体 偏 心率 较 大 或 轨道 倾角 较 高 时 ， 
特洛伊 ). 随 着 更 多 观测 的 开展 以 及 观测 技术 的 进 还 有 可 能 出 现 一 些 包含 以 上 两 种 轨道 类 型 的 混合 
步 ， 越 来 越 多 的 行星 特洛伊 被 发 现 . 截至 目前 , 已 轨道 . 
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Fig.1 The equilibrium points (L,—L;) and three types of co-orbital orbits (tadpole, horseshoe and quasi-satellite) in the planar 


circular restricted three-body problem 


一 最 早 人 们 将 处 于 太阳 -行星 的 暗 时 形 轨道 上 的 。 ”也 称 轨道 共振 ), 特洛伊 受到 来 自行 星 的 强烈 摄 动 
之 。 天体 称 为 特洛伊, 其 实在 一 定 扰动 时 山 峙 形 轨道 和 ”这 种 1:1 MMR 对 特洛伊 的 动力 学 演化 以 及 长 期 稳 
CO 马蹄 形 轨 道 会 相互 转变 , 现在 马蹄 形 轨道 上 的 小 天 定性 有 重要 影响 传统 的 Laplace 摄 动 函数 展开 中 
O 体 往往 也 称 为 特洛伊 . 特洛伊 具有 十 分 复杂 、 有 趣 ”包含 Laplace 系 数 ， 该 系数 在 两 个 天 体 与 中 央 天 体 
已 。 的 运动 特征 , 它们 的 来 源 和 去 向 也 可 能 成 为 研究 太 “距离 相同 时 是 发 艇 的 ,因而 不 能 处 理 11 MMR: 另 


一 阳 系 早期 历史 的 重要 线索 回 一 方面 ， 摄 动 函 数 展开 往往 基于 偏心 率 e 和 轨道 倾 
= 从 方法 来 看 对 特洛伊 和 共 轨 运动 的 研究 以 数 。。 角 ; 这 两 个 小 量 , 难以 处 理 e 和 i 较 大 的 情况 "2. 

CD。 值 模拟 居多 . 一 些 工作 针对 个 别 特洛伊, 对 观测 获 对 高 轨道 倾角 问题 , 可 以 针对 特定 的 倾角 构造 
CN 得 的 轨道 进行 数值 模拟 以 分 析 它 们 的 轨道 稳定 性 、 。 相应 的 摄 动 函数 tel， 例 如， 对 于 极 轨道 上 的 半 人 


Al 推测 其 来 源 , 例如 针对 新 发 现 的 地 球 特洛伊 名 和 海 马 座 小 天 体 和 海王 星 外 天 体 (Trans-Neptunian Ob- 
: 星 特洛伊 了 3 的 细致 分 析 . 另 一 些 工作 则 对 具体 jects), 它们 的 倾角 接近 90°, 可 以 将 e 和 cosi 作 为 小 
某 个 行星 的 三 角 平 动 点 附近 的 稳定 区 域 作 全 局 性 ERF. 对 于 道行 运动 ， 倾 角 接 近 180*， 则 将 e 和 
的 探讨 , 例如 最 引 人 关 注 的 木星 特洛伊 8 1 、 海 王 cos(i/2) 作 为 小 量 进行 展开 . 进而 , 对 于 任意 倾角 ， 
星 特洛伊 HL 五 电 和 地 球 特洛伊 03 均 的 情况 . MEK 可 以 设 定 一 参考 倾角 ii, 使 得 sin(i 一 i,) Av) at, 并 
批 太 阳 系 外 行星 系统 的 发 现 , 也 有 一 些 工作 不 局 限 以 此 进行 摄 动 函 数 展开 99. 利用 这 种 方法 构造 的 
于 有 具体 的 太阳 系 行星 , 而 是 考察 质量 参数 更 一 般 的 摄 动 函数 ,， Namouni 等 人 计算 了 2:1、3:1、1:2 等 共 
系统 中 三 角 平 动 点 的 长 期 稳定 性 9. 一 些 工作 更 进 振 的 共振 宽度 以 及 共振 中 心 的 位 置 23. 但 他 们 的 
一 步 , 不 局 限于 限制 性 三 体 问题 模型 , 而 是 研究 太 摄 动 函数 中 仍然 使 用 了 Laplace 系 数 ， 因 而 并 不 能 
阳 系 外 行星 系统 中 以 特洛伊 轨道 构 型 存在 的 行星 处 理 1:1 共 振 的 情况 . 为 了 处 理 共 轨 运动 问题 , 人 们 
的 运动 特征 19， 甚至 是 处 于 特洛伊 轨道 的 类 地 行 选择 不 使 用 Laplace 系 数 的 摄 动 函 数 展开 方式 ， 例 
Æ [17-18] 1 (它们 因 质 量 较 小 难以 直接 探测 ) 的 运动 及 guNamouni?4], Namouni##?5], Nesvorny“#?4); 


= i 


| a 


其 可 观测 信号 . 过 半分 析 的 方法 对 系统 的 哈密 顿 直接 进行 数值 平 
除了 数值 方法 之 外 , 也 有 不 少 工作 利用 简化 的 均 , 研究 e 和 i 较 大 时 共 轨 运动 的 长 期 演化 . 

限制 性 三 体 模型 ， 通 过 分 析 或 半分 析 模型 来 研究 上 述 工作 往往 聚焦 于 解决 特定 轨道 参数 范围 

特洛伊 的 动力 学 演 i 20), 在 太阳 -行星 -特洛伊 组 内 的 问题 , 而 较 少 对 三 维 连 续 参数 空间 中 的 共 轨 运 


成 的 系统 中 , 特洛伊 与 行星 的 平 运动 频率 相同 , 发 动 进行 系统 的 描述 . 本 文 利 用 新 近 出 现 的 一 种 摄 动 
生 1:1 平 运动 共振 (mean motion resonance, MMR, 函数 展开 方法 2 来 处 理 空间 圆 型 限制 性 三 体 模 型 
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中 1:1 MMR 问题 . 该 方法 不 使 用 Laplace 系 数 , 可 以 首先 将 该 摄 动 函 数 以 e 为 小 量 进行 展开 (为 保 
1 效 地 描述 共 轨 运动 ， 同时 该 展开 方法 适用 于 任 证 收敛 性 , 要 求 e < 0.6627), 然后 定义 

意 轨道 倾角 , 仅 对 偏心 率 有 一 定 限制 (小 于 0.6627， 

以 保证 展开 式 收敛 性 )， 本文 将 在 偏心 率 范围 e < Gae oee (4) 
(0,0.5), 倾角 范围 i € (0°,90°) 的 参数 空间 内 分 析 人 


k 轨 运动 , 特别 是 三 角 平 动 点 附近 的 运动 情况 . 本 
文 结构 如 下 : 第 2 节 简 单 介 绍 该 摄 动 函数 展开 方法 
以 及 本 文 所 用 的 数值 模型 ， 第 3 节 利 用 上 述 摄 动 函 
数 展开 以 及 数值 积分 计算 共振 中 心 随 e、i、 近 点 
角 距 w 等 轨道 根 数 的 变化 以 及 晴 蛙 形 轨道 的 发 生 范 
Hl; 最 后 在 第 4 节 总 结 本 文 结论 . 


2 方法 与 模型 
我 们 在 圆 型 限制 性 三 体 问题 模型 之 下 讨论 三 

角 平 动 点 附近 的 运动 ， 该 模型 中 两 个 主 天 体 ( 本 文 
中 分 别 为 太阳 和 木星 ) 在 万 有 引力 作用 下 沿 圆 轨道 
相互 绕 转 ， 而 零 质 量 的 第 三 体 ( 即 特洛伊 天 体 ) 在 
主 天 体 的 引力 场 中 运动 . 在 处 理 1:1 MMR 问题 的 
时 候 ， 我 们 希望 能 够 把 影响 小 天 体 运 动 的 摄 动 函 
P e,i,w,Q2,0) 的 函数 ， 其 中 a 为 小 天 体 
的 轨道 半 长 径 、9 为 升 交点 经 度 、o 为 1:1 MMR 的 
Lie FA, o 一 般 定义 为 小 天 体 与 行星 的 平 经 度 之 差 : 


= 


=] 


画 


g=A-N = (Wt 04 M)- (w F 01'4 M’), 


(1) 
星 的 


其 中 和 为 小 天 体 的 平 近 点 角 ， 
轨道 根 数 . 

在 限制 性 三 体 模型 中 , 小 天 体 围绕 第 一 主 天 体 
ee TA 离 二 体 
轨道 , 因为 引力 可 写成 摄 动 势 的 梯度 形式 , 所 以 可 
引入 摄 动 函 数 RR: 


带 撒 号 的 则 为 生 


1 
R=Gm (~~ eos) f (2) 
其 中 G 是 万 有 引力 常数 , r>、” 分 别 是 小 天 体 和 第 二 
主 天 体 m/ 到 第 一 主 天 体 m 的 距离 , % 是 小 天 体 和 7m/ 
相对 于 m 形 成 的 夹 角 , 而 人 是 小 天 体 和 mi/ 之 间 的 距 
离 , 表达 式 为 : 


A = (r? +17? — arr”? cosy) a (3) 
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其 中 a、a/' 分 别 是 小 天 体 和 mv 的 轨道 半 长 径 ,， a = 
a/a' 为 两 者 之 比值 . 显然 , O € [0,1], Ha = a 时, O 
= 1/2. EMO, = 2Q0/(1 + ao)2( 其 中 ao 为 初始 
半 长 径 之 比 ), 此 时 5 = 9 - 6. 是 一 个 小 量 , 可 将 摄 
动 函数 以 6 为 小 量 进行 展开 . 具体 展开 过 程 参见 文 
献 [27]， 本 文 仅 列 出 展开 式 对 快 变量 X 作 平均 之 后 
的 摄 动 函 数 形式 (为 简便 起 见 仍 记 作 RR): 


> > 


即 之 0 k20 


t l—t 


™~ ~ 
> >D Kox 
0 t1>0 


t220 


R= 


M 


m,(l—2tı—2t2) 
2t—l—2t,+2te 


Xo, 


pay Vi 
fav) 


sin” 


[Dn fala)] TUDE 
cos |(2t — l — 2t + 2te)o + 2(1 — t — 2t2)w] — 
cm (5) [Xi (e) cos? (é/2) cos(c) 十 
Xt (e) sin? (i/2) cos(o — 2w)] , 


(5) 


tpn, ky g bhom ty ty tol ARAM FP hy, 
和 jaax 分 别 为 e 和 6 的 展开 截断 阶 数 , 本文 分 别 取 
为 4 和 30，Xs"(e) 为 汉 森 系数 ，D 为 导数 算 子 ， 下 


标 n 代 表 对 a 的 n 阶 导数 , 且 
q 
ful) = Gy gy (6) 
而 ro 的 表达 式 如 下 : 


(—1)" m+- (2k — 1)! 


k\(q 
k)! C) 0 x 
n\/t\/l-t\/l @k-4 
m ti to t (1 = @,)1/2+k i 
其 中 () 为 二 项 式 系数 . 
可 采用 如 下 一 组 了 


Ko = 


(7) 


E 则 变量 来 表达 上 述 摄 动 函 


64 卷 天 
数 : 
P=Va, Q =o, 
P, = Ja(V1—e? -1), Qo =u, n 
P; = vVa(v1 - cosi- 1), Q3 =, 
P =X +P, Qs =, 
SOHN A 5 iT B-PAPEN SGC Ae. 进而 1:1 


MMR 系统 平均 化 之 后 的 哈密 顿 函数 (为 方便 起 见 ， 
以 常用 轨道 根 数 表 达 ) 为 : 


(9) 


在 该 平均 化 后 的 系统 中 , 不 显 含 Q3 和 Qa, 两 者 
为 循环 坐标 ， 与 之 共 斩 的 角 动 量 声 和 已 为 运动 积 
分 , 其 中 常数 已 给 出 了 系统 长 期 演化 过 程 中 a、e、 
;的 耦合 关系 . 本 文 则 在 研究 共振 时 标 内 三 角 平 动 
点 附近 的 运动 情况 , 而 w 的 变化 远 慢 于 共振 角 的 变 
化 , 因此 大 多 数 时 候 可 将 w 近 似 为 常数 , 哈密 顿 系统 
即 降 为 一 自由 度 的 可 积 系统 . 

除了 利用 上 述 摄 动 函数 展开 和 哈密 顿 形式 进 
行 半 分 析 的 计算 , 我 们 在 本 文中 还 利用 数值 方法 对 
运动 方程 进行 直接 积分 . 本 文 使 用 了 Swifter 数 值 积 
分 软件 ?中 的 Helio 算 法 5 积分 测试 粒子 的 轨道 . 积 
分 时 长 经 多 次 试验 选 为 50 倍 行星 的 轨道 周期 ， 积 4 
步 长 为 1/20 行 星 的 轨道 周期 , 积分 过 程 中 控制 精度 
为 测试 粒子 的 位 置 , 速度 相对 误差 小 于 10- 22, 中 心 
天 体 为 太阳 , 行星 在 圆 轨 道上 , 行星 的 半 长 径 为 单 


1 
H* ~ 94 — va — Ra, e,i,0,w). 


位 长 度 , 太阳 和 行星 总 质量 为 单位 质量 , 其 中 行星 


质量 为 10-3Me, 相当 于 木星 , Ms 为 太阳 质量 


共振 结构 
平面 圆 型 限制 性 三 体 模 型 中 ， 两 个 三 角 平 动 
点 DJ、Z75 即 为 1:1 MMR 的 两 个 共振 中 心 ， 对 应 的 
k 振 角 分 别 为 c = 60° Flo = 300°. Lyf Ls Fo = 
180°SE MT PK, ALLA SCT OL, PE AZ, 
Ls BAU AY ToL AY DAI RTT IAS 
3.1 ”共振 中 心 

平 动 点 对 应 着 摄 动 函 数 的 极 小 值 , 因 


2https://www.boulder.swri.edu/swifter/ 


3 


而 它们 的 


文 
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位置 可 以 由 摄 动 函数 展开 的 表达 式 求 出 . 实际 上 当 
小 天 体 的 轨道 不 是 共 面 的 圆 轨道 时 , 共振 中 心 的 位 
扁 离 60° 425]， 所 以 车 要 讨论 三 角 平 动 点 附近 
:1 MMR, 首先 需 确 定 共振 中 心 的 精确 位 置 . 
[前 所 述 , 共振 中 心 的 位 置 可 由 摄 动 函数 的 极 
小 值 位置 给 出 . 我 们 在 图 2 中 给 出 了 小 天 体 初 始 偏 
心率 eo = 0.1、 轨 道 倾角 io = 30?、 近 点 角 距 wo 
60°? 时 计算 共振 中 心 的 一 个 例子 . 图 中 摄 动 函数 
的 计算 值 ( 实 线 ) 随 初始 共振 角 oo 变 化 , 其 极 小 值 的 
位 置 以 虚线 标 出 ， 两 条 线 分 别 对 应 Za 和 Z5. 同时 
可 以 看 到 RR 的 变化 曲线 关于 co = 180? 对 称 , 在 o = 
180° 处 达到 局 部 极 大 值 , 该 处 对 应 着 平 动 点 Ls. 当 
给 定 小 天 体 的 (eo,io,wo) 时 , 我 们 可 以 通过 如 图 所 
示 的 方法 , 求 得 摄 动 函 数 极 小 值 的 位 置 , 该 位 置 即 
对 应 着 共振 中 心 . 


x10 


analytical 


numerical 


120. 180 300 


oo!” 
图 2 摄 动 函数 尽 ( 实 线 , 左 纵 轴 ) 以 及 半 长 径 振 幅 Aa( 散 点 , 右 纵 轴 ) 随 
* 振 角 co 的 变化 


240 


Fig.2 The perturbation function R (solid line, left ordinate) 
and semi-major axis amplitude Aa (scattered dots, right 


ordinate) versus the resonant angle oo 


另 一 方面 , 位 于 共振 的 小 天 体 其 共振 角 c 和 半 
长 径 a 发 生 秤 动 ， 而 在 共振 中 心 处 它们 的 秤 动 幅 度 
最 小 ， 因 而 我 们 也 可 以 用 数值 方法 模拟 小 天 体 的 
运动 ， 并 通过 检查 它们 的 秤 动 幅度 的 方法 找到 共 
振 中 心 . 图 2 中 的 散 点 在 对 数 坐 标 上 显示 了 我 们 用 
数值 方法 积分 运动 方程 所 得 到 的 半 长 径 的 秤 动 
(Aa). 同样 地 , Aa 在 oo = 60? 和 co = 300° HTE) 
到 极 小 , 与 计算 摄 动 函数 极 小 值 的 分 析 方 法 得 到 上 
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结果 一 致 . K, 在 io = 70?"-80? 附 近 达 到 极 值 (c。 ~ 56°). 小 天 


通过 上 述 两 种 方法 可 以 得 到 共振 中 心 (以 该 处 体 的 轨道 偏心 率 较 小 (eo。 = 0.1) 时 , 数值 结果 与 分 
URAC HIE NR, 以 下 记 作 Gc) 的 位 置 随 着 小 天 体 析 结 果 吻 合 得 很 好 ; 偏心 率 略 大 (eo = 0.3) 时 , 两 种 
轨道 根 数 eo、io、wo 等 的 变化 情况 , 我 们 总 结 了 几 方法 的 结果 略 有 差别 . 分 析 方 法 是 经 过 对 摄 动 函 数 
组 典型 根 数 情况 , 显示 于 图 3. 图 中 不 同 的 线 型 代表 的 平均 化 之 后 的 情况 , 而 数值 方法 没有 平均 化 的 过 
不 同 的 参数 下 数值 积分 的 结果 和 利用 摄 动 函 数 的 JE, 显示 了 有 瞬时 根 数 的 情况 . 受到 展开 阶 数 以 及 摄 
分 析 结 果 . 图 3 (a) 是 给 定 轨 道 偏 心率 情况 下 共振 中 动 函 数 收 敛 性 的 影响 , 偏心 率 较 大 时 两 种 方法 所 得 
心 随 倾角 ij 的 变化 情况 . 我 们 看 到 ce 先 减 小 后 略 增 的 结果 之 间 误 差 会 增 大 , 但 仍 在 合理 范围 之 内 . 


(a) (b) 


ey =0.1, numerical 


ig = 10°, numerical 


70 prann 


a e =0.3, numerical ----- ig =30°, numerical 
oo eo =0.1, analytical 80| --- io = 10°, analytical 
65 oN *  e9 =0.3, analytical °, analytical 
S 70 A 
CY 
C ‘ J 
A 60-——__ ee 
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(c) 1 四 
eo =0.1, io = 30°, numerical 70.5} ----- eo = 0.3, io = 10°, numerical 
iW 
_ 59.5| -一 - analytical * analytical 


图 3 共振 中 心 的 位 置 随 eo、io、wo 的 变化 


Fig.3 The variation of resonance center with respect to eo, io and wo 


图 3 (b) 显 示 了 给 定 轨道 倾角 io = 10" 和 30? 时 ， 给 定 偏心 率 和 轨道 倾角 的 情形 下 共振 中 心 随 近 点 
k 振 中 心 随 偏心 率 eo 的 变化 . 此 时 我 们 看 到 ce 随 偏 角 距 的 变化 . 在 eu = 0.1. io = 30°? 和 eo = 0.3、 
心率 的 增加 而 显著 地 变 大 , 在 较 小 的 偏心 率 (io。= io = 10? 两 种 情形 下 , 共振 中 心 o6 的 变化 范围 均 仅 
10°) 情 况 下 , 共振 中 心 可 达到 o。 ~ 85°. 分 析 方 法 与 有 1° 左 右 , 表明 近 点 角 距 值 wo 对 共振 中 心 位 置 的 影 
数值 方法 的 符合 情况 与 图 3 (a) 类 似 一 在 低 倾角 和 响 不 大 , 也 同时 表明 小 天 体 的 轨道 运动 几乎 不 受 近 
小 偏心 率 情 况 下 基本 完全 一 致 ， 随 着 i 和 eo 的 变 大 点 角 距 的 影响 . 所 以 为 简单 起 见 , 在 本 文 以 下 部 分 ， 
两 者 略 有 差别 , 在 我 们 关心 的 轨道 根 数 范围 内 , 这 若 无 特殊 情况 , 我 们 固定 一 个 wo 值 , 不 失 一 般 性 , 均 
种 差别 可 以 接受 . 取 为 60". Namouni 等 人 利用 Hil 三 体 问题 模型 给 出 

类 似 图 3 (a)-(b), 为 了 检查 初始 近 点 角 距 对 小 了 Ls4、Ls 点 位 置 的 解析 表达 式 B9, 发 现 共振 中 心 
天 体 运 动 的 影响 , 我 们 在 图 3 (c)-(d) 中 显示 了 两 组 的 位 置 主要 受 小 天 体 偏心 率 和 轨道 倾角 的 影响 , 而 


& 
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近 点 角 距 对 ce 影响 不 大 , 与 图 3 的 结果 一 致 . 

由 上 述 讨论 可 知 , 计算 摄 动 函数 的 极 小 值 即 可 
获得 共振 中 心 的 位 置 . 固定 近 点 角 距 为 60" 之 后 , 我 
们 在 偏心 率 -轨道 倾角 的 平面 (eo,io) 上 计算 了 共振 
中 心 的 位 置 , 并 展示 于 图 4. 图 4 中 空白 区 域 为 赂 时 
轨道 消失 的 区 域 , 而 虚线 为 数值 方法 给 出 的 量 蛙 
轨道 存在 的 边界 . 


图 4 (eo, io) PHE 


# 振 中 心 的 位 置 分 布 . KERRO. 一 60° 的 大 小 . 


Fig.4 The resonance center on (eo, io) plane. The grayscale 


indicates a, — 60°. 


在 图 4 中 , 我 们 可 以 看 到 在 轨道 倾角 较 小 时 , pE 
着 偏心 率 的 增加 共振 中 心 逐渐 远离 c。 = 60° 而 向 
La Ji ABB). 在 给 定 轨道 偏心 率 的 情形 下 , 共振 中 
心 的 偏离 又 随 着 轨道 倾角 的 提高 而 略 有 缩小 . 从 摄 
动 函数 的 角度 分 析 , 共振 中 心 的 这 种 偏离 都 是 使 得 
在 共振 中 心 处 摄 动 函数 值 尽量 地 小 . 而 从 几何 构 型 
上 看 , 这 样 的 偏离 使 得 小 天 体 在 椭圆 轨道 上 运动 时 
避免 过 于 接近 摄 动 天 体 (行星 ). 
3.2 ”共振 宽度 

当 偏 心率 和 轨道 倾角 为 零 时 , 共振 中 心 的 小 天 
体 处 于 受 力 平衡 状态 , 其 共振 角 和 半 长 径 取 为 平衡 
位 置 的 固定 值 . 当 小 天 体 的 初始 位 置 偏 离 平 衡 点 、 
或 初始 轨道 偏心 率 和 倾角 不 为 零 时 , 共振 角 和 与 之 
罗 的 轨道 半 长 径 即 以 一 定 的 振幅 振荡 . 人 们 将 平 
运动 共振 中 能 够 保持 共振 运动 构 型 的 天 体 半 长 径 
的 变化 范围 称 为 共振 宽度 , 共振 宽度 反映 了 共振 影 
响 的 范围 , 可 以 表征 共振 的 强度 . 与 半 长 径 振 幅 变 
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化 同时 发 生 的 还 有 共振 角 的 振幅 , 因此 后 者 有 时 也 
作为 共振 宽度 的 指标 , 也 可 以 通过 共振 角 的 行为 判 
断 运 动 与 共振 中 心 的 偏离 . 本 文 定义 小 天 体 运 动 中 
半 长 径 a 的 上 下 边界 (运动 的 边界 ) 为 共振 宽度 . 
3.2.1 共振 宽度 的 计算 方法 
我 们 已 经 展示 , 在 固定 了 近 点 角 距 w 之 后 ， 哈 
密 顿 函数 (9) 式 刻画 的 是 一 自由 度 的 系统 , 一 对 共 
斩 变 量 已 、Qi 构 成 该 哈密 顿 系 统 的 相 空间 . 在 共 
振 角 c 和 半 长 径 a 的 平面 内 作出 哈密 顿 函 数 等 能 线 ， 
就 可 以 了 解 该 共振 的 全 局 几何 性 质 . 以 eo = 0.1、 
io = 30°. wo = 60° 的 情形 为 例 , 我 们 作出 哈密 
顿 等 能 线 ， 如 图 5 所 示 . 图 中 实 线 代 表 哈 密 顿 等 能 
线 , 背景 灰 度 代 表 哈 密 顿 绝 对 值 的 大 小 , 颜色 越 深 
该 值 越 大 , 反之 则 越 小 . 虚线 标 出 了 两 种 运动 类 型 
分 界线 的 位 置 , 内 部 为 围绕 共振 中 心 的 秤 动 , 对 应 
HULA IS, 而 最 外 部 o 作 循环 运动 . 在 两 者 之 间 ， 
o 以 180° 为 中 心 秤 动 , 振幅 小 于 360°, 是 马蹄 形 轨道 
运动 . ape Pui AB LE, 马蹄 形 轨道 在 相 图 上 所 
占 的 区 域 非常 小 , 说 明 它们 发 生 的 几率 很 小 , 也 不 
是 本 文 研究 的 重点 , 因此 图 中 未 标 出 马蹄 形 轨道 的 
区 域 . 
图 5 中 哈密 顿 绝对 值 极 小 值 处 为 共振 中 心 的 位 
置 , 两 个 共振 中 心 分 别 在 和 Ls 点 附近 , 方形 点 标 
出 的 是 c = 60? 附 近 的 共振 中 心 , 点 线 标 出 的 即 为 
晴 星 形 轨道 的 共振 宽度 Aa, 其 计算 公式 如 下 [20, 271; 


Aa = asep — 00 = YÈ (R(x) — Ro]? , (10) 


上 式 中 asow 为 分 界线 上 对 应 着 共振 中 心 位 置 的 半 长 
径 值 , 它 满足 H*(asoy,o0s) = H*(a0,0u), 而 cs 和 ou 
分 别 代表 稳定 和 不 稳定 的 平衡 点 的 共振 角 ，R(o;) 
和 RR(o) 则 为 对 应 的 摄 动 函数 值 . 显然 , 此 处 稳定 平 
衡 点 即 为 前 节 求 出 的 共振 中 心 cs = ce， 而 不 稳定 
平衡 点 为 Ls 点 附近 的 平衡 点 ,对 应 图 2 中 摄 动 函 数 
在 oo = 180? 附 近 的 极 大 点 ,可 以 通过 与 计算 ce 类 
似 的 方法 求 得 . 利用 (9) 式 和 (10) 式 即 可 计算 小 天 体 
的 共振 宽度 Aa, 我 们 称 这 种 方法 为 分 析 方 法 . 男 外 ， 
小 天 体 的 共振 宽度 也 可 以 由 数值 方法 积分 试验 体 
轨道 运动 方程 并 寻找 共振 中 允许 的 最 大 半 长 径 偏 
离 的 方式 来 计算 , 我 们 称 这 样 的 计算 为 数值 方法 . 
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图 5 ”哈密 顿 等 能 线 图 (eo = 0.1, io = 30°. wo = 60°) 


Fig.5 Hamiltonian level curve for 
30°, wo = 60° 


eo = 0.1, to 


3.2.2 ”轨道 倾角 的 影响 
图 6 总 结 了 以 分 析 方 法 和 数值 方法 计算 共振 
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宽度 的 一 些 例子 ， 其 中 初始 偏心 率 eo 分 别 为 0.1、 
0.3 和 0.5, 我 们 显示 了 共振 宽度 随 轨道 倾角 的 变化 . 

图 6 中 横 坐 标 为 试验 体 的 初始 半 长 径 co， 纵 坐 
标 为 粒子 的 初始 倾角 io. 试验 体 共振 角 取 为 共振 中 
心 的 位 置 co = oo, 由 前 一 节 的 计算 结果 给 出 ， 而 
近 点 角 距 固定 为 60°. 根据 (10) 式 计算 得 到 不 同 i 于 
的 Aa, 在 图 中 点 出 1 土 Aa, 由 红色 点 表示 . 在 用 数值 
方法 计算 共振 宽度 时 , 半 长 径 在 (0.9,1.1) 范 围 内 均 
匀 取 50 个 点 , 而 轨道 倾角 在 (0°, 90°) 范 围 内 取 90 个 
点 , 用 数值 方法 模拟 这 些 轨道 至 50 倍 轨道 周期 , 并 
记录 轨道 运动 过 程 中 半 长 径 的 变化 范围 (图 6 中 以 
颜色 显示 , 由 深蓝 到 浅 蓝 表示 半 长 径 的 变化 范围 由 
小 到 大 ). 由 图 6 容易 看 到 , 在 共振 中 心 处 这 些 试验 
体 的 半 长 径 振幅 很 小 , 离 中心 处 越 远 振幅 越 大 , 并 
且 在 特定 的 位 置 振幅 发 生 一 个 跳 变 (运动 由 师 蛙 形 
轨道 变 成 马蹄 形 轨道 或 其 他 类 型 轨道 , 图 中 由 最 右 
侧 的 黄色 点 标 出 , 这 些 点 通过 将 马蹄 形 轨道 的 平衡 
位 置 代入 (10) 式 计算 得 到 ). 该 位 置 就 是 星 蛙 形 轨 道 
振幅 最 大 处 , 其 与 共振 中 心 的 距离 也 就 是 共振 的 宽 
度 . 


eo=0.5 


1.05 


图 6 不 同 偏 心率 下 共振 宽度 随 轨道 倾角 的 变化 


Fig.6 The relationship between resonance width and inclination for several eccentricities 
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图 6 显示 数值 方法 给 出 的 共振 边界 与 由 分 析 方 
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心率 变化 的 情况 总 结 于 图 7 中 . 在 不 同 轨道 类 型 之 


法 给 出 的 边界 (共振 宽度 ) 非 常 好 地 相互 吻合 . 在 偏 
心率 较 小 时 (图 6 中 eo = 0.1、0.3), 分 析 结 果 与 数值 
结果 吻合 较 好 ; 而 当 eo = 0.5 时 , 分 析 方 法 与 数值 方 
法 的 吻合 情况 稍 差 , 主要 是 因为 偏心 率 较 大 时 分 析 
方法 使 用 的 摄 动 函数 收敛 性 变 差 . 同时 , 值得 注意 
的 是 , 相 比 于 小 偏心 率 的 情形 , 在 eo = O.5 EN MRE TZ 
轨道 出 现 的 区 域 大 幅 地 减 小 , 在 轨道 倾角 30° 以 上 


间 ， 共 振 和 角 振 幅 Ao 会 发 生 跳 变 , 可 以 比 半 长 径 振 
幅 Aa 更 好 地 区 分 轨道 类 型 ， 因此 为 了 更 清晰 地 刻 
画 边界 , 我 们 在 数值 结果 中 以 Ac 标识 和 区 分 不 同 
的 轨道 类 型 . 另外 值得 注意 的 是 , 考虑 到 高 偏心 率 
时 摄 动 函数 展开 式 收敛 性 的 限制 ,图 7 中 分 析 结 果 
只 显示 到 eo = 0.5, 而 数值 结果 则 到 eo = 0.99 为 止 . 
如 图 7 所 示 ， 随 着 偏心 率 增加 ， 曙 蛙 形 轨道 的 


的 区 域 , 星 昱 形 轨道 基本 消失 . 整体 上 看 , 倾角 对 共 
振 宽 度 的 影响 不 大 , 蜂 蚜 形 轨道 的 共振 宽度 随 倾角 
增加 而 略微 收 罕 . 
3.2.3 ”轨道 偏心 率 的 影响 

与 轨道 倾角 不 同 , 小 天 体 的 初始 偏心 率 显 著 地 
影响 1:1 共 振 的 稳定 性 与 轨道 构 型 89, 偏心 率 越 大， 
Lie He PE eS. 与 图 6 类 似 , 我 们 选择 了 典型 的 三 
个 轨道 倾角 io = 0*、30°* 和 60°, 将 共振 宽度 随 偏 
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图 7 不 同 轨道 倾 


a 


范围 逐渐 减 小 直至 在 特定 偏心 率 处 完全 消失 . 在 平 
面 模型 中 (io = 0°), WEH en = 0.8844 
WA, 被 围绕 Ls 的 李 雅 普 诺 夫 轨 道 ([Lyapunov or- 
bit) 取 代 . HAM TELE A ZB, 在 此 类 轨道 类 型 
的 外 侧 较 罕 的 区 域内 , 则 有 马蹄 形 轨道 (图 中 以 星 
形 表示 ) 出 现 . 另外 如 前 所 述 , 我 们 看 到 分 析 方 法 的 
结果 只 到 eo = 0.50 为 止 , 且 分 析 结 果 与 数值 结果 给 
出 的 共振 宽度 基本 一 致 . 
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共振 宽度 随 偏 心率 的 变化 . 从 左 至 右 分 别 对 应 io = 0°, 30°, 60° 的 情况 . 图 中 星 形 (*)、 叉 号 (X) 和 加 号 (十 ) 分 别 代 表 马 蹄 形 轨 


i. Labi Lyapunov ili. wR ES HE TEKRARI, 红色 点 为 分 析 方法 给 出 的 共振 宽度 , 灰色 虚线 为 晴 星 形 轨道 消失 时 的 偏心 率 . 背景 


颜色 为 数值 方法 计算 获得 的 共振 


振幅 人 Ao. 


Fig.7 The resonance width versus eccentricity for io = 0°,30° and 60°. Stars (x), crosses (x) and pluses (+) denote horseshoe 


orbits, the Lyapunov orbit around Ls and the combination of tadpole and quasi-satellite orbits!?°], The red points mark the 


resonance width obtained by analytical method, while the grey dotted line marks the eccentricity when tadpole orbits vanish. 


Color in background indicates the libration amplitude of resonant angle (Ac). 
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在 倾角 较 高 (io = 30°. 、60°) 时 , HRMS AZ Eit 
会 在 更 低 偏心 率 ( 分 别 为 eo = 0.49、0.44) 处 消失 . 
但 与 平面 模型 不 同 的 是 ， 较 高 轨道 倾角 的 情况 下 ， 
T gC Z Ab RS FZ Se AA H E E eS FF 
准 卫 星 形 混合 而 成 的 轨道 类 型 (图 中 加 号 )， 而 马蹄 
形 轨道 仅 出 现在 eu < 0.10 时 晴 蛙 形 轨道 区 域外 侧 
的 狭 窗 范 围 内 . 图 中 虚线 上 方 区 域 , 大 偏心 率 区 域 
出 现 的 小 振幅 (Ac) 运 动 是 准 卫 星 轨道 . 另外 , 轨道 
倾角 较 高 时 蜂 虹 形 轨道 仅 出 现在 较 低 偏心 率 处 , 此 
时 分 析 方 法 所 揭示 的 赂 时 形 轨道 共振 宽度 与 数值 
方法 得 到 的 结果 基本 一 致 . 

鉴于 分 析 方 法 在 较 低 偏心 率 的 情况 下 获得 的 
晴 蛙 形 轨 道 的 共振 宽度 与 数值 方法 给 出 来 的 结果 
基本 一 致 , 我 们 利用 分 析 方 法 计算 (eo,zo) 平 面 上 是 
星 形 轨道 的 共振 宽度 , 并 将 结果 总 结 在 图 8 中 , 此 图 
即 给 出 了 蜂 虹 形 轨道 存在 的 范围 . 图 中 颜色 代表 共 
振 宽 度 (Aa) 大 小 , 而 空白 处 不 存在 晴 蛙 形 轨道 

由 图 8 可 知 ， 在 较 高 轨道 倾角 或 较 高 偏心 率 
的 情况 下 不 存在 晴 时 形 轨道 . 实际 上 图 6 中 eo = 
0.50. io > 30° 的 区 域 以 及 图 7 中 io。 = 30*、eo > 
0.49 和 io = 60°, eo > 0.44 的 区 域 , 小 天 体 仍 在 1:1 
MMR, 其 运动 也 仍然 是 稳定 的 (Aa 或 Ao 的 规则 分 
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意味 着 该 处 的 轨道 仍然 具有 规则 而 稳定 的 特征 )， 
只 是 轨道 的 类 型 变 成 了 准 卫星 轨道 . 本 文 主要 关注 
晴 蛙 形 轨道 的 情况 , 但 为 了 说 明 运 动 类 型 的 变化 ， 
我 们 在 图 9 展示 了 初始 偏心 率 eo = 0.50, 而 倾角 分 
别 为 io = 20" 和 40?* 时 的 哈密 顿 相 图 . 


0.0500 
0.0475 
0.0450 
£ 
0.04252 
z 


0.0400 


\ | 


i 
| 
| 


0.2 0.3 
eo 


0.0375 


0.0350 


图 8 (eo, io) 平 面 上 共振 宽度 的 分 布 


Fig.8 The resonance width in (eo, tio) plane 


io = 40°, œo = 60° 
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图 9 初始 倾角 20?( 左 ) 和 40? ( 右 ) 的 哈密 顿 等 能 线 图 , 初始 偏心 率 为 eo = 0.5 


Fig.9 Hamiltonian level curves for io = 20° (left) and 40° (right), with eo = 0.5 for both 
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图 9 中 的 哈密 顿 等 能 线 图 明显 不 同 于 小 轨道 偏 。 ” 统 ) 轨 迹 在 (o, a) 平 面 上 的 投影 , 即 为 图 中 橙色 点 , 而 
ne ih 一 0.1 时 的 图 5) 小 偏心 率 情形 ”黑色 点 为 对 应 测试 粒子 的 初始 ele 可 以 看 到 
F, RA La Ls XTM KIS E ZA LXT h E H 测试 粒子 的 运动 轨 迹 基 本 符合 哈密 顿 等 能 线 的 走 
不 动 点 较为 显著 , 它们 分 别 是 映 蚜 形 轨道 和 马蹄 形 ” ”向 , 说 明 该 哈密 顿 系 统 仍 可 a 
轨道 的 共振 中 心 ; 而 在 较 大 偏心 率 eo = 0.50 的 情 。 ”运动 . 值得 注意 的 是 , 图 9 中 左 侧 (io = 20*) 显 示 的 
BR, o = 0? 附 近 新 出 现 的 不 动 点 控制 了 相当 大 测试 粒子 的 轨迹 同时 包围 了 晴 糙 形 轨道 和 准 卫 星 
的 区 域 , 该 不 动 点 即 为 第 三 主 天 体 ( 行 星 ) 附 近 运 行 ”轨道 的 共振 中 心 , 这 正 是 图 7 中 以 “+4” 标示 的 “混合 
的 准 卫 星 轨道 的 共振 中 心 . 同时 , 对 于 给 定 的 偏心 ”轨道 ”. 
率 , 哈密 顿 相 图 还 随 着 轨道 倾角 而 变化 . 图 9 中 io = 
20* 时 ,我 们 仍然 可 以 看 到 原 斑 、 疙 附近 的 晴 时 形 «3-3 近 点 角 距 的 影响 
轨道 的 共振 中 心 , 而 当 i0 = 40° 时 , 原 到 附近 的 共振 图 9 已 显示 疡 和 7 附近 的 共振 结构 不 对 称 , 元 
中 心 消 失 了 , 由 此 处 出 发 的 轨道 将 围绕 c ~ 0 的 共 附近 的 晴 旦 形 轨 道 在 更 小 的 轨道 倾角 条 件 下 先 消 
振 中 心 秤 动 , 成 为 准 卫星 轨道 . 此 外 , Spit Rempel RDN. 发 生 这 一 现象 的 原因 主要 是 我 们 计算 时 给 
轨道 的 共振 宽度 类 似 , 准 卫星 的 共振 宽度 也 可 以 通 。 ” 定 了 wo = 60° 这 一 不 对 称 的 值 (wo = OOM FL. L 
过 分 析 方 法 计算 , 我 们 在 图 6 中 出 现 准 卫星 轨道 的 。” 才 是 对 称 的 ). 虽然 图 3 表明 偏心 率 较 小 时 晴 蛙 形 轨 
子 图 (eo = 0.50) 上 已 经 相应 标 出 了 分 析 方法 计算 得 。” 道 的 共振 中 心 随 wo 的 变化 较 小 , 但 在 偏心 率 较 大 
到 的 准 卫星 轨道 的 发 生 范 围 (图 中 橙色 点 组 成 的 边 “时 ,共振 结构 会 随 着 wo 的 取 值 不 同 而 发 生 改 变 29, 
Fi). 本 文 不 对 此 作 深 入 研究 , 只 在 图 10 中 作出 eo = 0.1. 
我 们 在 图 9 中 二 加 了 积分 测试 粒子 运动 方程 而 0.5 两 种 偏心 率 条 件 下 wo 分 别 取 0*、60° 和 120° 的 哈 
得 到 的 真实 运动 (而 非 截断 且 平 均 化 后 的 哈密 顿 系 BWA, 以 显示 wo 的 取 值 对 轨道 构 型 的 影响 . 
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振 结 构 的 影响 


Fig.10 The influence of wo on resonance structure 
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由 图 10 中 左 侧 子 图 可 知 , wo = 0° 时 , 不 管 偏心 距 会 引导 小 天 体 进入 不 同 的 共振 轨道 构 型 . 
率 大 小 ,Za 和 万 附近 的 运动 区 域 都 是 相互 对 称 的 . 以 上 结论 , 都 分 别 利用 摄 动 函数 展开 的 分 析 方 
而 图 10 中 上 一 行 子 图 则 表明 , 偏心 率 较 小 (eo = 0.1) 法 和 直接 积分 运动 方程 的 数值 方法 给 出 , 两 者 相互 
时 ，wo 的 不 同 取 值 对 相 空 间 结 构 的 影响 几乎 可 以 对 比 验 证 , 也 证 实 了 该 摄 动 函数 展开 方法 处 理 共 轨 
忽略 ，L4 和 Ls 附近 的 运动 仅 有 微小 区 别 ; 而 在 较 运动 的 适用 性 . 
大 偏心 率 (图 10 中 下 一 行 子 图 , eo = 0.5) 的 情况 下 ， 
wo = 60°? 时 (下 中 图 ) 7 附近 的 晴 昱 形 轨道 区 域 消 
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On the Co-orbital Motion of Any Inclination 


LU Rui! LEI Han-lun'”? ZHOU Li-yong!? 


(1 School of Astronomy and Space Science, Nanjing University, Nanjing 210023) 
(2 Key Laboratory of Modern Astronomy and Astrophysics, Ministry of Education, Nanjing 210023) 


ABsTRACT The co-orbital motion occurs when a celestial body (e.g. an asteroid) shares the same semi- 
major axis with the perturbing object (e.g. a planet), and thus they are in a 1:1 mean motion resonance. 
Trojan asteroids in the tadpole orbits around planets in the Solar System are these co-orbital objects. 
The motion and origin of some Trojan asteroids, particularly those on high-inclination orbits, are still not 
fully understood. In this paper, a newly developed perturbation function expansion, which is applicable to 
the 1:1 resonance, is used to investigate the co-orbital motion in three-dimensional space. The position of 
the resonance center and the resonance width are calculated for different initial orbital elements, and the 
relationship between the orbital types and the initial orbital elements is analyzed. The results obtained 
by the analytical method are compared with and verified by the results from numerical simulations. A 
panorama of the co-orbital motion in the wide orbital elements space is obtained. 


Key words celestial mechanics: circular restricted three-body problem, methods: analytical and numer- 
ical, co-orbit 
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